Superharmonic Josephson relation at 0-/7r-junction transition 
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Critical current was recently measured near the transition from to 7r-contact in superconduc- 
tor/ferromagnet /superconductor Josephson junctions. Contrary to expectations, it does not vanish 
at the transition point. It shows instead a tiny, though finite, minimum. The observation of 
fractional Shapiro steps reenforces the idea that the vanishing of the main sinusoidal term in the 
Josephson relation gives room to the next harmonics. Within quasiclassical approach we calculate 
the Josephson relation taking into account magnetic scattering. We find that the observed minimum 
is compatible with the value of the second harmonics expected from the theory. 

PACS numbers: 74.50.+r 



According to textbooks equilibrium supercurrent, /, in 
a tunnel-barrier Josephson junction depends sinusoidally 
on the phase difference, <p, between the superconducting 
leads: I = I c sin <f>, where I c > is the so-called criti- 
cal current of the junction. In his seminal work of late 
seventies Bulaevskii et al predicted Q that the sign of 
I c can change (or, equivalcntly, a shift of ir appear in 
the argument of the sine) in the presence of magnetic 
impurities within the tunnel barrier. Soon after, Buzdin 
et al. 0, Q suggested that such a junction, convention- 
ally called now the 7r-junction, can be realized in a hy- 
brid structure where the tunnel barrier is replaced with 
the ferromagnetic metal. While predicted theoretically, 
experimental realizations of 7r-junctions remained long 
unobserved. Indeed, superconductivity and magnetism 
compete; thus conventional ferromagnets would strongly 
suppress supercurrent. The first recent successful realiza- 
tion of a 7r-junction 0,0] utilized the so-called weak fer- 
romagnets, and the observation of a non-monotonic de- 
pendence of I c as a function of the temperature 0, Q and 
on the thickness of the ferromagnetic layer [{| served as 
the first evidences of the actual realization of 7r-junctions. 
Moreover in the former case, the existence of the temper- 
ature T*, where I c reached a minimum with the vanish- 
ing magnitude, allowed for a precise identification of the 
transition point. 

These spectacular observations of 7r-junction behaviors 
remarkably confirmed original predictions of Q, 0, 0] and 
yet posed new puzzling questions. The first was that 
the observed amplitude of the current in the junctions 
appeared two orders of magnitude smaller than that ex- 
pected from the theory. H. Sellier et al 0] proposed that 
magnetic impurities in the ferromagnet could be the ori- 
gin of this effect, and indeed it was recently shown for the 
sin (p component Q that magnetic impurities can lead to 
a noticeable reduction of the critical current if one as- 
sumes somewhat artificial uniaxial distribution of mag- 
netic disorder. It remains however to understand the 



effect of a more realistic disorder distribution. 

Another puzzle concerns the form of the phase-current 
relation at the transition point T* . Quite generally, the 
phase relation has to be periodic in </>. This does not rule 
out the possibility of the second or even higher harmon- 
ics: I = I\ sin <j) + li sin 20 + . . . , which indeed appear 
in Josephson junctions formed by point contacts, con- 
strictions, or non-equilibrium normal metals 0,0]. Only 
in the non-equilibrium case was the I2 component ob- 
served [f2T ] . However the amplitude of higher order com- 
ponents in the magnetic Josephson junctions was long 
considered too low for being observed. Note now that 
at T* the coefficient of the first harmonics (Ji) vanishes 
and, therefore, the higher harmonics become dominant. 
In Ref. [lflj the measured critical current (the maximum 
of the absolute value of the current-phase relation) does 
not vanish at T*, but passes through a minimum. This 
fact, together with the observation of fractional Shapiro 
steps, indicates that the observed current is in fact the 
I2 component. A phenomenological decoherence time has 
been proposed in Ref. 0] to fit both I\ and I2 ampli- 
tudes at zero temperature. On the other hand no I2 
component was detected in Ref. 0] ■ 

In this Letter we develop a microscopic theory enabling 
the quantitative derivation of the full current-phase rela- 
tion in the regimes corresponding to actual experiments 
of Refs. 0, 0, 0| . We solve the resulting equations 
numerically without restriction on the values of the pa- 
rameters. Extracting the exchange field and magnetic 
scattering time from the published data on the temper- 
ature dependence of J c , we estimate the expected mag- 
nitude of 12- We find that the predicted values agree 
favorably with those observed in Ref. , while the ex- 
pected magnitude of I2 for the sample of Ref. 0] is 
too small to be observable. Note that thickness inho- 
mogeneities in the ferromagnetic layer could give an al- 
ternative explanation [lfij ] of the experimental result of 
Sellier et al. [lflj . A powerful and microscopic approach 



2 



to superconductivity in disordered metals is offered by 
the quasiclassical theory in a form described in [rR Il8j . 
The theory can also describe ferromagnetism, by inclu- 
sion of an exchange field acting on conduction electrons. 
This was done for instance in Refs. 0, |2(|, where the 
spin-orbit coupling with impurities was also included. 
However, for the weak ferromagnet Cu^Nii-^ used in 
experiments of 0, 0, 0| , the spin-orbit coupling is 
expected to play a minor role. The more important ef- 
fect should come from the strong inhomogeneities of the 
magnetic field on both the microscopic- (magnetic im- 
purities) and mcsoscopic scales (randomly oriented mag- 
netic domains). Our theory takes this effect into account. 

The metallic ferromagnet is described by the following 
Hamiltonian: 



M+ U ) S SS f - h.a s 

2m 



where V's( r ') and ipl( r ) are annihilation and creation op- 
erators for electrons having spin projection s along the 
z direction, m is the effective electron mass, and /_t is 
the Fermi energy (h = 1). The disorder potential U(r) 
describes the interaction of electrons with nonmagnetic 
impurities and is characterized by the correlation func- 
tion: U(r)U(r') = 5(i — r') j (2-kvt), where r is the clastic 
mean free time and v is the density of state at the Fermi 
level per spin. The upper bar stands for disorder averag- 
ing. The exchange field h(r) acting on the electron spins 
may originate, for instance, from contact interaction be- 
tween conduction electrons and localized impurity spins. 
We do not consider the question of the microscopic ori- 
gin of h(r), but restrict ourselves to setting its statistical 
properties only. Namely, we take its average to be spa- 
tially uniform: h(r) = hz, with h proportional to the 
magnetization of the ferromagnet. The fluctuating part 
is characterized by correlation functions: 
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2-kvtSL 



5 a0 S(r-r'), (1) 



for a,/3 = x,y,z. Here characterizes mean free time 
due to magnetic impurities. In the following, we also 
assume rotational symmetry around z, thus rf n = t* . 

In order to describe the proximity effect in the fer- 
romagnet, it is convenient to introduce thermal Green's 
functions £„ s ,„' 8 '(r,r',-r) = -(T T ^ ns (r,T)¥ n , s ,(r',0)) 
in the Nambu(n)-spm(s) space, where = ip s and 
^2s = V'-s- The equation of motion for the Matsubara- 
transformcd disorder-averaged Green's function, Q, is de- 
rived from the Hamiltonian and reads: 



2m 



/' 



ha. 



Q = \. (2) 



Here, ui n are Matsubara frequencies at temperature T. 
The non-magnetic disorder leads to the usual forme of 



the self energy Ei (r,w„) = (2-kvt) 1 T z Q(r 7 r 7 uj n )T z . The 
magnetic disorder gives instead: 

S 2 (r,cj„) (27n/T£)-V*S a S(r,r,u; n )S a T>. (3) 

a— x,y ,z 

Here S = (a x , a y , <j z t z ), a a and r a are Pauli matrices in 
spin and Nambu spaces, respectively. 

Now we define the quasiclassical Green's function 
g(r,uj n ) = [i/TTv)T z Q(r,r,uj n ) which, in the diffusive 
limit, obeys the equation 



-DV(gVg) + [u> n r z - ihr z a z 



y— 



S a gS a ,g] = 0, 

(4) 

with D the diffusion coefficient, and the normalization 
condition g 2 = 1. Symmetry properties of the Hamilto- 
nian further constrain the form of g. Specifically: (i-a) 
By the invariancc under the rotation around the average 
magnetization axis i, we find that g is block diagonal in 
spin space; we thus define the two matrices in Nambu 
space, g+ and <?_, as the two non- vanishing upper and 
lower components, respectively, (i-b) From the invari- 
ance under the rotation over an angle 7r around the x 
(or y) axis and simultaneously change of sign of h, we 
find g+(h) = T z g-(—h)r z . (ii) By time-reversal symme- 
try, we find g±(r,u n ) = -T z g±(r, -w„)V z . Finally the 
16 correlation functions introduced in Eq. (??) are not 
independent, since the representation is redundant. This 
gives: (iv) g + (r,w n ) = -T x g_(r,u) n )*T x . 

Exploiting these properties, we parameterize the com- 
plete Green's function in terms of g+ 1 which in its turn 
is completely determined by complex functions, 6 and m 



.9+ = 



cos 9 sin Oe" 1 
sin 0e~ lTI — cos 6 



(5) 



Then, Eq. © yields: 



£>V(sin 2 9Vr]) + — sin6>sin6>* sin (77 - 77*) (6a) 
-DV 2 6 + Dcos9 sin 6( V77) 2 



+2(u> n — ih) sin^ -I sin cos ( 

T z 



sin cos ( 



cos 6* sin 0* cos(?7 — 77*)] (6b) 



These equations constitute the main analytical result of 
our work. Note that in non ferromagnetic superconduc- 
tors, symmetry properties (i-b) and (iv) for h = imply 
that 9 and 77 are real. Then Eqs. © only depend on the 
effective magnetic scattering time l/r m = \jrf n + 2/r^, 
in agreement with Abrikosov-Gor'kov theory for mag- 
netic impurities [2ll ]. By contrast, in ferromagnetic su- 
perconductors, magnetic disorder can be characterized 
by two scattering times: rf, = tK. and r* [13, ■ In 
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model 


order 


h (meV) 


l//lT m 


p{%) 


Rh (nV) 




1 


4.6 


9.3 


25 


1.08 




2 


49 


1.1 


8 


0.25 


(isotropic) 


3 


83 


0.3 


10 


0.07 




1 


15 


2.7 


12 


0.10 


T"m = OO 


2 


61 


0.8 





0.16 


(uniaxial) 


3 


86 


0.2 


8 


0.08 



TABLE I: Parameters of the fit (shown in Fig. for the 
isotropic case). We took for the sample A = 1.3 meV, D/L 2 = 
1.13 meV for L = 17 nm, T* = l.l K. p = 07/(207 + j b L) is 
the ratio of the barrier resistance to the total resistance R of 
the junction in its normal state. 



Ref. |jj the uniaxial disorder was considered. In our no- 
tation this corresponds to ^ = = 00, = r m . This 
hypothesis simplifies greatly the solution of Eq. © , since 
9 and 77 are no more coupled to 9* and rj* . The physical 
reason for this simplification is that magnetic scattering 
does not couple the spin up and spin down populations. 
However it seems more realistic that the magnetic disor- 
der is also able to flip the spin of conduction electrons. In 
this sense the opposite limit is to consider a completely 
isotropic disorder: = = r m = 3r m . In the fol- 
lowing we thus concentrate on this case. We shall also 
discuss briefly the uniaxial one for comparison. 

In order to determine the Josephson relation, we as- 
sume that the ferromagnet is a layer of the length L along 
the i-axis. We thus need to solve Eqs. JSJl with appro- 
priate boundary conditions [3| at x = ±L/2: 



07 yjujl + A 2 2 



(7a) 



V0 



^_ lb uj n sin 9 - A cos 9 cos(?7 T f) ^ 



°7 



\A4 + A2 



Here, 07 is the conductivity of the ferromagnet, and jb 
is the barrier resistance per unit area at the contacts 
(taken to be identical for simplicity) between the ferro- 
magnet and the superconducting leads. We also assume 
that temperature dependence of the superconducting gap 
in the leads, A, is simply given by conventional BCS the- 
ory. The supercurrent is given by 



2nGfLT \ - _ r . 9 n 
/ = - 2^ Re [sin 2 9\7i]] , 

W„>0 



(8) 



where G/ is the conductance of the ferromagnetic metal. 

We use now the above equations for fitting the exper- 
imental data of Ref. 0, ^| . For a given set of param- 
eters, namely the Thouless energy D/L 2 , A(T = 0), T, 
76, h, T m , we obtain the current-phase relation by solv- 
ing numerically the system of differential equations © 
and calculating the current through Eq. |(SJ|. One can 
then extract the first two harmonics, I\ and Ii- (Higher 
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Temperature (K) 

FIG. 1: Fit to experimental data for the critical current 
of Ref. [l3j| (boxes), for the three solutions reported in the 
table of the isotropic model. The three solutions gives nearly 
indistinguishable curves. Inset: temperature dependence near 
T* of the calculated I c — max^,(|/(0)|), Ji (dashed), and I2 
(dotted), for the first solution. The minimum of I c coincides 
with 12- 



harmonics near T* are much smaller than I2). We begin 
with the data of Ref. 01, Fig. 2, concerning a sample 
with a ferromagnetic layer of L = 17 nm. The length L, 
the superconducting gap, and the temperature are known 
experimentally. The interface resistance is more difficult 
to measure. The authors of Ref. 0] give an estimate of 
30% of the total resistance of the junction in its normal 
state, R 0. 

For a given value of jb one can find the pairs of values 
(h, t) that satisfy the two equations Ii(h,T m ,T*) = 
and Ii(h,T mi Ti) = 1^ , where T\ is a temperature dif- 
ferent from T* and I^ xp is the corresponding experimen- 
tal value for I±. We find that only three pairs of val- 
ues satisfy this constraint. We order them by increasing 
value of h. One can show that the n-th solution refers to 
a junction where, as a function of the length of the sam- 
ple, other n — 1 zeros are predicted for L < 17 nm. We 
optimize then this first estimate of the parameters by in- 
cluding 7t, as a fitting parameter for the full experimental 
curve. The solutions are given in Table [I] For compari- 
son fitting parameters for uniaxial magnetic disorder are 
also given. 

In all cases we have a good fit to data (see Fig. [IJ. 
Thus, the quality of the fit is not a sufficient criterium 
to discriminate between the three possibilities (for each 
model). One argument in favor of the first solution (for 
isotropic model of magnetic disorder) is the agreement of 
the fitting parameter jb with the estimated value in the 
experiment. A second one is the dependence of T* on 
the length that increases with the order of the solution. 
The range of L for which < T* < T c is about 1 nm for 
the first solution, and about 0.4 — 0.3 nm for the second 
and third one. The first case compares better with the 
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FIG. 2: Length dependence at T = 4.2 K of I c (full line) and 
I2 (dashed line) for the parameters obtained with the fit with 
the isotropic model. We show the solution 1 (upper panel) 
and 2 (lower panel). The zero at L = 17 nm is visible in both 
cases. Solution 2 displays a second zero for L < 17 nm. 



a function of the temperature |3j ■ 

In conclusion, we have presented a development of 
the quasiclassical theory of superconductivity taking into 
account magnetic scattering in the presence of an ex- 
change field. We have used our model to extract the 
exchange field and scattering times from the tempera- 
ture dependence of the critical current in superconduc- 
tor/fcrromagnet /superconductor junctions. With these 
parameters we have calculated the second harmonics at 
the vanishing value of the first component which agree 
favorably with the experimental findings. 

Note added: after the completion of this work we be- 
came aware of related work by A. Buzdin |2,"| where I2 
is calculated near the critical temperature. 

We are pleased to thank A. Buzdin for illuminating 
and important discussions. This work was supported in 
part by the U.S. Department of Energy office of Science 
under the contract W-31-109-Eng-38. 



experiment [l3j, where incertitude on L is about 1 nm, 
and 7r-contacts were observed for L = 17 — 19 nm [||. 

We consider now the second component, 1%. 

In Ref. ^3 the minimum value of RI2 is 0.5 nV, and it 
falls between the first and the second predicted value for 
isotropic model (cf. Tab. [IJ. In both cases we thus find 
that the amplitude of the second harmonics is compatible 
with the observed one. We also find a strong temperature 
dependence of I2; the values presented in the table at 
T* = 1.1 K, arc reduced by a factor 10 at 5 K. This 
dependence can explain the much smaller value for RI2 
observed in the 19 nm sample of Ref. [l3T | . In comparison, 
the uniaxial model gives a much smaller value for RI2 ~ 
0.1 nV for all three solutions. 

We repeated the fitting procedure on the data of Frolov 
et al. [15J, where no second harmonics is observed at 
T* . We found again that the data can be compatible 
with either a second or a first zero, but in both cases 
RI2 < 10~ 10 mV, thus below the observation threshold. 

We finally discuss the length dependence of the first 
and second harmonics for the fitted values of the pa- 
rameters (see Fig. |2J. As anticipated, Ii(L) displays an 
oscillating behavior. One can clearly sec in Fig. 0that I\ 
for solution 1 and 2 vanishes once and twice, respectively, 
for L < 17 nm (solution 3 is not shown). A more unex- 
pected result is the oscillatory behavior for h{L), that 
shows a remarkable doubled periodicity with respect to 
Ji(i): Between two zeros of the first harmonics we al- 
ways observed two zeros of the second harmonics. This 
means that the sign of I2 remains always positive when 
I\ vanishes. Therefore we find that the transition from 
0- to 7r-contact is always discontinuous. We cannot rule 
out, however, that I2 may be negative at T* in some 
other region of the parameters space. This would imply 
that the transition from to 7r-contact is continuous as 
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